ASYMPTOTIC SHAPES WITH FREE BOUNDARIES 



PHILIPPE DI FRANCESCO AND NICOLAI RESHETIKHIN 

Abstract. We study limit shapes for dimer models on domains of the hexagonal lattice 
with free boundary conditions. This is equivalent to the large deviation phenomenon for 
a random stepped surface over domains fixed only at part of the boundary. 



Introduction 

In this paper we study limit shapes for random tilings by rhombi with free boundary 
conditions. 

Recall that there is a local bijection between 3D-partitions, tilings of a plane (of a 
triangular lattice on a plane) by rhombi, and dimer configurations on the dual hexagonal 
lattice |6j. Another important local bijection is between lattice paths and rhombus tilings. 
It is illustrated on Fig{TJ 

The limit shape phenomenon is studied well in dimer models on hexagonal lattices and 
for corresponding tiling models in the case when the boundary of the domain consists of 
long pieces parallel to the axes of the triangular lattice. There are also some results for 
the case when the boundary is a fixed "zig-zag" path. But in all cases the shape of the 
boundary is fixed. 

The goal of this paper is to study the limit shapes when part of the boundary is a random 
variable. We call this "free boundary" . An example of a free boundary for a cut hexagon 
region is shown on Fig. [TJ The end points of lattice paths on the Noth-East boundary are 
not fixed, which corresponds to variable configurations of boundary tiles. 

In this paper we assume that configurations (tilings, 3D-partitions, paths) are distributed 
uniformly, or weighted with the probability 

(0.1) Prob(S) oc g 151 

where q is a real positive number, and \S\ is the volume under the height function corre- 
sponding to the configuration. 

The computations are based on the Gessel-Vionnet formula for the number of lattice 
paths [3]. In the last section we compare our results with the analysis of limit shapes of 
height functions for dimer models based on the Burgers-Hopf equations [I]. 

In Section [1] we compute the boundary of the limit shape for cut hexagons (see FigJT]) 
of large size when the lattice paths are distributed as in (10. II) . We assume that as the size 
of the cut hexagon increases, q is approaching 1 as q = exp(e) where - is proportional to 
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(0,n) 




Figure 1. A typical tiling of the cut hexagon nx k x n, with free boundary 
conditions on the cut. We have represented the tiling on the square lattice. 
One keeps track of the two non-square tiles by drawing a segment in the 
middle. The segments form non-intersecting lattice paths, starting at points 
(i — 1, 1 — i), i = 1,2, ...,k and ending on the cut x + y = n, at points 
(rrii, n — m;), < rrii <n + k — 1. 

the linear size of the system. In Section [2] we extend these results to the case when there 
are "frozen" triangles along the free boundary. In particular, when two frozen triangles 
are next to the ends of the free boundary, the problem includes the non-cut hexagon as a 
particular case Figj5l In Section [3] we compare our results for the full hexagon with the 
paper [2] where they were derived originally. Finally Section H] contains the comparison 
of our results with the derivation of the limit shape from the Hopf-Burgers equation [I]. 
In Section [5] we use a symmetry principle and the results obtained in previous sections to 
obtain the limit shape for the TSSCPP problem [TJ. 

The work of N.R. was supported by Danish National Research Foundation through the 
Niels Bohr initiative, by the NSF grant DMS-0601912, and by DARPA. P.D.F. is supported 
in part by the ANR Grant GranMa, the ENIGMA research training network MRTN-CT- 
2004-5652, and the ESF program MISGAM. 

1. Case of a cut hexagon 

In this section we will focus on random tiling with weights flO.ip . The goal is to find the 
asymptotic distribution of boundary positions of paths in the continuum limit i.e. when 
k, n — > oo and q = exp(±e) — > 1 such that ke — > f3 and ne — > a are finite. 

We will not give a rigorous proof of the existence of the limit shape, i.e. that the 
random boundary partition converges in probability to the limit shape function. We will 
give reasonable standard arguments, based on the variational principle, that it does, and 
compute the limit shape explicitly. 
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Figure 2. The area under the path (green boxes) plus that in the top tri- 
angle (pink boxes) accounts for the partition function Q (l) q for lattice 
paths from (0, 0) to (5, a — b), with a weight q per box in the corresponding 
slice of plane partition. 



1.1. g-enumeration. The generating function for the number of tilings viewed as plane 
partitions, with a weight q per box (with green top) can be expressed as: 

( m j- i + i )( m j - i ) ( n 

(1.1) Z k , n (m\q) = det q 



l<i,j<k \ m j — i + 1 

where the g-binomial reads (£) = ^ jjjj^gj and the g-factorial [a] q \ = Yl^ =1 (l — 

This formula is explained in Figf2J The g-binomial (fy enumerates the lattice paths 

from (0, 0) to (b, a — b) with a weight q A , where A is the (discrete) area under the path, 
namely the number of unit squares in the domain delimited by the axes x = 0, y = a — b 
and the path. This gives the correct weight q per box in the plane partition, within the 
slice containing the path, provided one adds up the contribution of the boxes in the domain 
delimited by x = 0, y = a — b and x + y = a (pink boxes), with area (*) , hence the additional 

pref actor gv 2 ) . 

The formula (II. ip factors out a Vandermonde type determinant and it has the following 
multiplicative form: 

(1.2) Z M (m|g) =g m^^ 1 ^(^-^D A(0 f[ ,f n + fc ~^ ! 



J- [mi] q \[n + k-rrii-l] 
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where A(g m ) = rii<i J -<fc(? m< " Q mj )- 

The generating function is symmetric with respect to the mapping q — > q^ 1 

(1-3) Z k , n (m\q) = q s ^Z Kn (m\q) 

where rhi = n + k — 1 — m.j, Sk, n = ~ nk ( n + k ) _|_ _|_ ^ _ \y 

1.2. The continuum limit. 

1.2.1. Now we consider {m} random variable with 

Prob({m}) oc Z fc)n (m|g) 

and we will study this random variable in the limit when q = e ±£ , e — *> +0, n,k—>- oo, and 
a = ne, [3 = ke are fixed. 

Let us start with the case q = e~ £ . The following results can be easily extended to q = e e 
using the symmetry ( 11. 31) . 

When e — > and /x = me is kept finite we have: 

Log[m],! = -/ Log(l - e~ x )dx + 0(1) 
£ Jo 

Assume that as e — > 0, the density 

fe 

/i(t) = mie5(t — is) 

i=l 

converges to the density corresponding to the continuous function //(£) on the interval [0, 0\. 
Then, the leading term of the asymptotic of the generating function (11.21) is given by 

Z M (m|«) = e~^ +0 ^ 

where 

rP i rP f M*) /-A+l-Kt) 

= / -/i(*)(M(*)-2/9)+ / / Log(l-e- x )^+ / Log(l - e"*) cfe 

Jo ^ J Jq Jq 



P rP 

Log | e -M(s) _ 6 -mW |d s dt + C(a, P) 



1 

2 jo Jo 
where C(a,/3) reads: 

-/3 ra+p-X 



C(a, /3) = -- J dx J Log(l - e^) dy 



Assume that the function fj,{t) is different iable and consider the density p(/i): 

H\t) 
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In terms of p(p) the functional S reads as 



S[p] = J (^p(p-2(3) + L(p) + L(a + (3-p))p(p)dp 

j ras+P rc 

Vo Jo 



a+/3 r a+f3 

In |e _M — e~ u \p{p)p{v) dp dv 



2 

where L(x) = JTlr^l — e~*)<i£. 

In addition, since Y^!i=i 1 = we have 

(1.4) /?= / dt= / p(A*)d/x 



The functional S[p] is convex and therefore has unique minimum. 
It is easy to find the variation of this functional: 



ra+P 

[1.5) 5S[p] 



^p(p - 2(3) + L(p) + L(a + (3- p) 

»a+/3 r a+/3 



Sp(p)dp 

pa+p pa+p 

f f In \e-^ - e- au \p(u)du5p(ii)dfi 
Jo Jo 



The second variation is positive definite: 
6 2 S[ P ] 



pa+p pa+p 

+ + In |e" M - e~ v \5p{v)5p(p)dpdv 
Jo Jo 



Therefore the minimum of S[p] is achieved at its critical point. Let us find the critical point 
of S which satisfies the constraint (11. jj) . The convexity of S guarantees its uniqueness. 
The condition (jl.4p implies the constraint on the variation of p: 

-a+P 

Sp(p) dp = 



Taking it into account, we obtain the integral equation for critical points of S[p\: 

1 ra+p 

-p(p - 2(3) + L(p) + L(a + (3- p) - f In |e _M - e~ v \p{u) dv = const 

2 Jo 

We want to solve it on the space of functions p satisfying the constraint (11.41) . 
Differentiating this integral equation and taking into account (1 1.4ft we arrive at 

ln( ) = [ a+P piy) d(e-») 

After changing variables to u — e~ v and L = e _Q_/3 the equation for the density cr(u) 
p(u) becomes 

(1.6) H Z -^1) = I 1 ^± 

l-Z It z-u 
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where the function a satisfies the condition 

(i.7) r^ =/3 



u 

1.2.2. Now let us solve ( 11. 61) on the surface of the constraint ( 11.71) explicitly. 
Proposition 1.1. The unique solution to U.6]) subject to the constraint fll. 7| ) is 



9 / /sinh sinh (jB ~ m) 

(1.8) p(/i) = a(e^) = -Arctan e^J , 2 A . g 2 

7r \ V smh # smh % 



2 2 



where 



-b -2±a / ,«+^ / • ,0 • , 2a + /3 
e = e 2 cosh \ / smh — smh 



2 V 2 2 



- i _-2+£ f « + , /_.__,_ (3 _ z __,_ 2a + /3 



^1.9) e = e 2 ( cosh h y sinh — sinh 

2 V 2 2 

Proof. Assume the support of a is the interval [a, 6] C [L, 1] and consider the function 

/*o(«)du 

on the complex plane with the branch cut along [a, b] . It is a meromorphic function whose 
boundary values at the interval [a, b] satisfy the relation 

2 ln(— -) = F(z + iO) + F(z - zO) 
1 — 2 

and i* 1 ^) = as |z| — > oo. 

To solve this equation explicitly consider the function G(z) = , 1 =F(z). It is a 

meromorphic function on the complex plane with the branch cut along [a, b]. But now at 
the branch cut we have 

G {z + iO) - G(z - iO) = 2 \ = ln(^— -) 

and G(z) = 0{\) as |z| — > oo. From here we obtain: 

. . 1 /7 w 7T Z" 6 1 U — L du 

(1.10) F(z) = -y/(z-a)(z-b) Ti = =Log- 



7T 

As U| — ► oo this function behave as 



„ y/(u - a) (b - u) l-uz-u 



F W = - / 7? wt r Log^ ^ + O - 

- a) (6 - u) 1 - u z 
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This gives the first equation for a, b: 

1 u — L 
— ^^^=^^^= Log du = 

y/(u-a)(b-u) l-« 

It is easy to see that the integral vanishes when a + b = 1 + L. 

Changing the integration variable in ( 11.101) u = (b — a)s + a and taking into account the 
relation between a and b we obtain the following integral formula for F(z): 

1 n 77 7T f 1 1 T s + M ds 

F(z) = -y/(z-a)(z-b) / -== 



7T 

«— a r^^A A/f _ a— £ 



Log 



o v/s(l - s) & l + M-s £- 



where t = f-A and M 



6-o' l+L-2a' 

This integral can be computed explicitly using 

/ (w-wVtf + wUMl -«)*"" + w (. 



+ 



+ 1) ^/w{w - 1) 
and we obtain: 

. - , z — L . (z — a)(b — z) 

F(z) = In + MArctan J j J - 

1 — z y (a — L){1 — a) 

The second equation determining a and b is the normalization (jl.4p . It gives F(0) = —f3 
and therefore: 



ab 2 
This gives a and b: 



e 2 cosh \ / sinn — smh 



2 V 2 2 



Q + /3 



(1.11) 6 = e 2" I cosh 01 — — + A / sinh — sinh Q — ^ 

! 2 V 2 2 

The density p(/i) can be easily derived from F{z) as p(p) = imF(e~ M ) when ^ G [5, A], 
and 5 = —Log b, A = —Log a: 



;i.l2) p(/i) = a(e^) = -Arctan e" M 



2 ( sinh - „ ; sinh 



^l, (B-M) 



2 2 



7r \ V sinh ^ sinh -| 



□ 



Because the solution p described above is a local minimum of the convex functional S. 
it is also a global minimum. 

The density p(p) is plotted on Figj3] for A = 1/2,1,2,5 and (3 = 0,1/2,1,2,5. Here 
A = a/ (3. 
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X =2 =5 

Figure 3. Distribution of endpoints for A = = 1/2,1,2,5 and (3 = 
0, 1/2, 1, 2, 5 (from right to left, in each plot). 

1.2.3. The limit when a = f3X and a — > while A is fixed correspond to the limit shape 
for uniformly weighted lattice paths, i.e to q — 1. 
In this limit A = j3A , B = j3B where 



;i-13) B (X) = I±*_</±(1 + ^) 



(1-14) 4,(A) = ^ + V^ 1 + ^) 

The density of endpoints on the free boundary converges to 



here t = fj,/(3. 

1.2.4. Another interesting limit is when a — > oo while /3 is fixed. In this case p is given 
by the same formula as before but for the boundary points we have: 

1 



;i - V 7 ! - e-P) 



A - ^(i + v / T^) 
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Figure 4. The domain to be tiled is a half-hexagon nx k xn with the fixed 
interval of m square tiles along the boundary. 



Note finally that all the above is valid for a < as well. The transformation a 
amounts to reflecting the plane partition w.r.t. the axis y = x as in (11.31) . 



—a 



2. Case of a cut hexagon with forbidden intervals along the free 

boundary 

In the previous section we studied boundary partition for the distribution g' s L From 
now on we will consider the uniform distribution. All results can be easily generalized to 
q 7^ 1 but we will not do it here. 



2.1. Enumeration. We consider the domain of FigJH made of a cut hexagon n x k x n 
with a fixed interval of square tiles along the boundary. The effect of these tiles is that 
lattice paths are forbidden to have endpoints along a segment of the form {(i, n — i), i = 
Hi, Tii + 1, . . . rii + Tn — 1}, while ni + m + ri2 = n + k. 

It is clear that the counting formula is the same Gessel-Viennot formula as in the previous 
problem with q = 1: 



(2.1) 



Z k ,n( m ) 



det 

l<i,j<fc 



n 



7Tlj 



It can be written in the multiplicative form: 
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2.2 Z k , n (m)= \l {m i -m j TT 1 l - -y 

-"r- 1 - J/ ±± [mi\l[n + k — — 1J! 

The difference with the previous case is that positions of end points of lattice paths m« 
now can only be in [0, n + k — 1] \ [ni, ni + m — 1] = [0, ni — 1] U [ni + m, n + A; — 1]. The 
interval [ni, n\ + m — 1] is forbidden. 

2.2. Large size asymptotics. Here we will study the distribution of endpoints in the 
limit when k —>■ oo and n = Xk, n\ = X\k, n\ + m = = fak, for some finite A, Ai, 

and A 2 . 

The large k asymptotics of the number of lattice paths with endpoints fixed outside of 
the forbidden intervals is given by the action function 



(2.3) S x [p] = / (-/i(/i-2A) + L (/i) + L (A + l-/i))p(/i)d/i 

(2.4) — -+ 4 In \n — u\p{n)p(v) dfx dv 

2Jo Jo 



where L (x) = xLog(x) — x. 

This is a convex functional. The same arguments as in the previous section shows that 
its extremum is the unique minimum. Assume that the forbidden segment [Ai, A2] is in a 
generic position which means that the region where the density is not constant (not or 1) 
consists of two segments T = [a, b] U [c, d], such that 0<a<6<Ai<A2<c<c?<A + l. 
The density is constant outside of T: 

p(ji) = 0, pe[0,a]U[d,\ + l], p(/i) = l, n e [b, Ai] U [A 2 , c) 

The density p(p) which minimizes S[p] satisfies the integral equation 

g Log u{b-u){u-\ 2 ) = I p(p)dp 

(A + 1 — w)(Ai — u)(u — c) J T u — p 

u G T with the additional constraint: 

p(u)du + (Ai - b) + (c - A 2 ) = 1 

T 

The solution to the integral equation (12.51) with the additional constraints can be obtained 
in the similar way as before by converting the integral equation to the problem of finding 
the holomorphic function with given boundary values. For the resolvent 



F(z) 



p(u)du 



T 



z — u 
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we obtain the following integral expression 

(2.6) {z) = -\/{z - a)(z - b)(z - c){z - d) 

71 

(2.7) x / - du — T™- «(&"«)(« -A 2 ) 



y/(b - u)(u - a)(c -u)(d-u) z-u (A + 1 - u){\\ - u)(u - c) 

The end points a,b,c,d are determined by the asymptotic behavior F(z) = (1 + A 2 — 
Ai + b — c)/z + 0(1/ z 2 ) at large z. Indeed, if there were no conditions on a,b,c,d the 
resolvent would behave as z 2 for large z. Vanishing coefficients in z 2 , z, z°, and fixing the 
coefficient in z~ l give four equations. If the parameters A, Ai,A2 are in generic position, 
these equations have a real solution, which gives the desired values of a, b, c, d as functions 
of A, Ai, A 2 . 

The density is given by the imaginary part of the boundary value for the resolvent on T: 

(2.8) p(z) = Im F(z + iO) /vr, z G T 

For sufficiently small values of Ai the interval [a, b] may shrink to a point. The same can 
happen with the interval [c, d] for sufficiently large values of A 2 . Both intervals shrink to 
a point when A 2 — Ai = A. It is impossible to have A 2 — Ai > A because in this case the 
number of lattice paths incoming to the region would be greater then the number of paths 
outgoing from it. 

2.3. Case of many forbidden intervals. This is a straightforward generalization of the 
previous case to a succession of forbidden intervals along the free boundary. 

The number of lattice paths is still given by the formula (12. ip . The only difference is 
that all rrii should be outside of the forbidden intervals. 

Denote by [a*, the forbidden intervals i = 1, . . . ,1. The density of the outgoing paths 
is supported on U^[/3j_i, a*] where /3 = and c^ + i = A + 1. Let us say «j, $ are generic 
when the minimizer of the rate functional S[p] is constant on the following intervals: 

p{u) = l,ue U{t£[&_i, a*] U [bi-u ai_i], p{u) =0,n£ [0, ai] U [b l+1 ,\+ 1] 
The integral equation for the minimizer is 

Lo ( - irr — - hi ^^ ~ l \ = I dw p ^ 

»A + l-/i-f^ (aj_i - p)(cii - p) I J T p-w 



i=2 



where T = Uj[«i, 

The minimizer satisfies the constraint (the number of incoming paths is equal to the 
number of outgoing paths): 



z+i i 



i=2 i=l 
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The solution to this integral equation is given by the the resolvent function 

F(z)= 1 - fduT\^MkzM T og — TT(^-^)^-A-i) 

tt J T x . x y/(bi - u)(u - Oi) z - u A + 1-m- l .- l (ai_i - u) {ai - u) 

as usual: 

p(z) = lmF(z + iO)/ir, zeT 
2.4. The case of forbidden and fully packed intervals. 

2.4.1. Forbidden intervals are just one way to partially enforce boundary conditions along 
segments of the free boundary, where one imposes p = 0. Another special case is that of 
fully packed intervals, along which every point is an end point of a lattice path, hence p = 1 
is imposed. 

We cannot enforce a given type of tiles along a fully packed interval. However the 
geometry of tiles is such that each fully packed interval has a sequence of non-square tiles 
which start with horizontally tilted tiles at the top and then horizontal tiles turn at some 
point into vertical non-square tiles. The position of the turning point is random: it is not 
fixed by the requirement that the interval is fully packed. 

The counting formula for the number of lattice paths is the same as before and we can 
compute the limit density of end points of lattice paths along the free part of the boundary 
similarly to the case when we have only forbidden intervals. 

2.4.2. The problem of finding the limit density for one fully packed interval is very similar 
to the problem for one forbidden interval. Lattice paths come to the fully packed interval 
horizontally at the upper part of the interval and vertically at the lower part of the interval. 
Such a configuration induces forbidden intervals at both ends of the fully packed interval. 
Let l\ and Z 2 be end-points of the fully packed interval. In the limit when k, n,li,l2 — » oo 
such that X — n/k, Xi — li/k are finite the resulting density of lattice paths has the following 
structure 

p(p) = 0, p e [0,a] U [d, A + 1], p(p) = 0, fi e [6, Ai] U [A 2 , c] 

Here 0<a<6<Ai<A 2 <c<rf<l + A, p(p) = 1 when p 6 [Ai, A 2 ] and p(p) is a 
smooth functions when p G T = [a, b] U [c, d] . It satisfies the constraint 

p(p)dp + (A 2 - Ai) = 1 

The limit density on T is determined by the variational problem similar to the one described 
above and the solution is given by the resolvent function eq. (12.61) . The density is equal to 
the imaginary part of the resolvent on T eq. (12.81) . The endpoints a, b, c, d are determined 
by the condition F{z) = (1 + Ai - X 2 )/z + 0(1/ z 1 ). 

When a fully packed interval is near any boundary of the whole interval [0, A + 1], the 
lattice paths form a "fully packed region" along the corresponding side of the cut hexagon. 
In this case the problem of finding the limit density of paths is equivalent to the problem 
for a cut hexagon of smaller size. 
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2.4.3. The generic case of several fully packed and forbidden intervals corresponds to an 
alternance of the intervals, apart from one-another. All other situations are degenerations 
of this one. 

The equation for the density in this case is similar to the case when we have several 
forbidden intervals. Let [aii, @i\ be forbidden intervals and [A*, Vj\ be fully packed intervals. 
Assume that $ < Aj for all % — 1, . . . , n. When the boundaries of intervals are in generic 
position the density will vary smoothly in intervals [a*, bi], [c$, di], i = 1, . . . , n where • • • < 
bi < di < fa < Ci < di < \i < [ii < a i+ i < b i+ i < . . . , and d n < A + 1. The density is 1 
in fully packed intervals and in the forbidden intervals. In addition in is 1 in intervals 
[bi, aii], [Pi, Cj], and it is in the intervals [di, Aj], a i+ i], [0, ai], and [d n , A + 1]. 

The density function satisfies the constraint: 



n „ n 

/ p(^W + Y((ui-bi) + 
1=1 ^ T 1=1 



(ci - pi) + (vi - A^) 



where T = U™ =1 [aj,6j] U [ci,dj] The solution to the integral equation for the minimizer is 
given by the resolvent: 

F ( z j = 1 [ duf[ V( z - a i)( z - b i)( z - c i)( z - d i) 



X 



i \/{bi -u){u- a,i)(u - Ci){u - di) 
1 u -A- (u - bi)(u - Pi){u - Aj) 



z — u A + 1 — u (u — ai)(u — Ci)(u — //j) 



The constraint on the density is equivalent to the following condition on the asymptotics 
of F(z): 

n 

F{z) = (1 - ^((oi - bi) + (q - pi) + (ui - \)))lz + 0(11 z 2 ) 

which gives equations defining bi, q, dj. The density is given by the imaginary part of 
FO) on T as in eq. fl2T8|) . 



2.5. Case of two forbidden intervals at the corners. Here we consider the case of two 
forbidden intervals that touch respectively the top and the bottom of the free boundary. We 
will see that in this case the limit density can be computed explicitly in terms elementary 
functions. 

Assume that the upper interval has length m and that the lower one is of length p. Now 
the values mi G [0, m] and rrii G [n + k — p — l,n + k — 1] are forbidden, and hence the 
allowed range for endpoints of lattice paths is [m + 1, n + k — p — 2]. 

Now we will take k — > oo such that n = Xk, m = vk, n + k — p = 9k, which finite A, v, 9 
and will compute the minimizer of the rate function S[p\. 
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2.5.1. For generic is, 9 the minimizer of the rate functional is not constant on the interval 
[a, b] with is < a < b < 9, and 

p(u) = l,u e [is, a] U [b, 9] 
It satisfies the integral equation 

(2.9) Log ( ^- <>)(«-£ ) = /' dw £^L 

V(A + 1 -fi)(ji- is)(b-p)J J a p-w 

and the constraint 

(2.10) I p(w)dw = 1 - (a - is) - (6 - b) 

J a 

In this case the integrals can be computed explicitly. 
Proposition 2.1. The density is given by 



, n , s , s 2/. b — z la — is K lb — z la 

(2.11) p(z) = 1 + - ArctanW «/- Arctan \ J- 

7T\ V z — o V b — is V z — ay b 



z — a 19 — b I z — a / A + 1 — b 

+Arctan <i / \ / Arctan 



b — z V 9 — a \ b — z \ X + 1 — a 

(2.12) 

with 

(2.13) a = ^(U(\, is, 9) - y/V(X,u,9)) b = -^(U(\, is, 9) + y/V(\,u,9)) 

where U and V are the following polynomials: 

U{X,is,9) = (is + 9)(l - (is - 9) 2 ^) - (X + - is - 6)(1 + 2X + is + 9) +4is9 

V(X,u,9) = (l-is-9)(l-u + 9)(l + is-9) 

x(l + is-9 + 2A)(1 -is + 9 + 2A)(1 -U-9 + 2X) 

Proof. The resolvent for (12.91) is given by the integral 

_ VX™ ~ !) f 1 dv 1 Log f (P + v)v(t] - v) 



71 



o y/v(l -v) w-v \(<y-v)(5 + v)(l-v) 



where w = f— - and 

o—a 

n a X + l — a a — is 9 — a 

b — a b — a b — a b — a 

The normalization condition (12.101) for the density is equivalent to large z asymptotic 
behavior for the resolvent F(z) = 1+fe ~ a+I/-6> + 0(\). Vanishing of the coefficient in z and 
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the constant term of this asymptotic amounts to: 

(2.15) = I' - dV Log ( ^ + VHV - V) 



y/v{l ~ v) \(J-V)(5 + V)(1-V) 

b — a f 1 dv ( (3 + v)v(r] — v) 

(2.16) 1- (5 + r])(b-a) = / v Log ' 



Using identities 

(2.17) t * Log(/3 + u)du = 2vrLog f^+^SH) 

Jo ^u(l -u) V 2 J 

(2.18) jT 1 ^_A__ M Log(/? + u)du = \{[^~ V^+l) 2 + 2Log + ^F+l ^ ^ 
we obtain equations 

i = ^|2(5 + r / + i) + ( v ^TT- ^y + (^i-^y 

-(Vs + i - V6) 2 - ( v^FT - v^) 2 } 

The integral defining the resolvent can be evaluated explicitly using the following idenity 



Jw(w - 1) f 1 dv 1 A -\- v 
— / = Log — 

7T Jo a/ v(l - v) w - V B + v 



/~,^ t A + W /. /1-tU / 5 . 1-W / A 

(2.19) Lo Sd + 2« Arctan W \ - --Arctan 



B + w \ V w \ll + B V w \l + A 

It results in the formula: 



F(z) = Log ( ; ^ ± W ) WiV ~ W) , )+2i Arctan - A ~ w ■ / 5 



(7 — u>)(5 + w)(l — w) J \ V w V 1 + S 



1 — w (3 / u> / 77 — 1 / /7 — 1 

—Arctan -\ / * / + Arctan * / * / Arctan 



w y 1 + (3 \ 1 — w \ r] y 1 — u; y 7 

Using this expression for the resolvent, vanishing conditions (12.15p(12.16p can be written 



as: 



(2.20) a = -^(U(X, u, 9) - y/V(\,u,9)) b = -^(U(\, u, 6) + \/V{X, u, 6)) 
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where U and V are the following polynomials: 

U(X,v,9) = {v + 9)(l-(v-9) 2 ^ -(X + l)((l-v-9)(l + 2X + v + 9)+Av9 

v(\,u,e) = (i-u-e)(i-u + e)(i + u-e) 

x(l + v-9 + 2A)(1 - v + 9 + 2A)(1 - v- 9 + 2A) 
From the formula for the resolvent we can compute the density z G [a, b\. 



(2.21) p(z) = lmF(z + iO)/n = 1 + - [ Arctan \ - — - Arctan \ - — 

7r \ \ z — a \ b — v V z — a V b 



z-a 9-b jz - a / A + l 

+ Arctan -i / \ / Arctan 



b — z V 9 — a \ b — 2 V A + l- a 

(2.22) 

valid for u, 9 satisfying v < a < b < 9. Note that p(a) = p(b) = 1. □ 

2.5.2. When the size of any of the forbidden intervals in sufficiently small, it merges with 
the nearest end of the interval [a, b]. For example when 



9 C = Ul + X + u+ v/3(l + 2A) + (1 + A - 2z/) 2 ) 



6 ■ 3 

the lower forbidden interval connects to [a, b\. b = 9 C . When 9 < 9 C they remain connected, 
i.e. in this case we still have b = 9. 

The same holds for the upper forbidden interval, in which case when 



1 -(l + X + 9- v/3(l + 2A) + (l + A-20) 2 



c 3 

the upper forbidden interval connects with [a, b], a — v c and they stay connected for v < v c . 

The assumption v < a < b < 9 which we used in the previous section holds when 9 < 9 C 
and v > v c . 

In 9 > 9 C and v > v c the region where p{u) is not constant becomes [a, 9 C ] with 

b(l + X + u-b) 2 
a= P 

The density in this case can be obtained from (12.111) by setting b,9 —> 9 C . The result is 



™x /x 2/ 9 r — z a — v 

(2.23) p{z) = - Arctan 



7r \ \ z — a y 9, 



9 c -z fa j z-a /X + 1-9 C 

Arctan \ / , Arctan 



z — a y 9 C \ 9 r — z \ X + 1 — a 

It is easy to see that now the density satisfies p(a) = 1, p(b) = 0. 
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^.-x e 




e-x 



0<x< 




x-A. ! 



e+A-x 




Figure 5. We decompose the hexagon 1 x 9 x A into two cut hexagons with 
one side equal to 1, and with possible triangular holes at corners. Three cases 
must be considered. 

The formula above is valid as long as v > u' c , namely where 



v c \0=6 C 



1 + A - V1 + 2A 



When v < v' c the upper forbidden interval connects to the interval [a, b\. In this case 



1 + A- V1 + 2A 



1 + A + V1 + 2A 



2 ' 2 

The density inside of [a, b] can be obtained from ( 12. 23ft by taking a, v 



is' The result is 



%-z){z-v' c ) 



which agrees with (11.151) . 



Remark 2.2. The transition from free to trapped paths is similar to that found by Douglas 
and Kazakov in QCD. 

2.5.3. When 9 + v — 1 the interval [a, b] collapses into a point. 

3. Case of the full hexagon 

In this section we will reproduce the results of [2] for the limit shape of the hexagon 
tiling and will prepare the ground for the next section where we will use the symmetries 
of regions to find the limit density at free boundaries for more complicated regions then in 
previous sections. 

Cut the hexagon as it is illustrated in Figj5] into two cut hexagons. In this section we 
will study the limit density of tiles by studying the the limit density of lattice path points 
intersecting an arbitrary cut (dashed line) which is parallel to one side of the hexagon . We 
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assume that the hexagon has size k x n x p, and that the cut is at a distance m from the 
edge of length k. Our goal to pass to the limit k —>■ oo with n/k = A, p/k = 9, m/k = x 
and to compute the limit density of lattice paths through the cut. We will assume A > 9 
which clearly does not restrict generality. 

The cut separates the hexagon into two regions. For the purpose of counting lattice 
paths intersecting the cut we distinguish three possible cases of geometries of these regions: 

• the lower region is a cut hexagon, while the upper region is a cut hexagon with two 
forbidden intervals, 

• both lower and upper regions have one forbidden interval, 

• the lower region is a cut hexagon with two forbidden intervals, while the upper 
region is a cut hexagon. 

All three cases are illustrated on Figl5] and have been treated in the previous sections. 
The new ingredient is that we must now glue the two halves, by identifying the endpoints 
rrii <-> m' fc+1 _ i of both families of paths. 

In the large k limit the density of paths through the cut develops the limit shape which 
minimizes the action functional determined by the asymptotic of the counting formula or 
the lattice paths. The formula is the same for both pieces of the hexagon. The resulting 
action functional is the sum of two pieces. 

To describe the resulting action let us parameterize the endpoints of paths of the lower 
half of the hexagon by ^ = rrii/k . In all three cases of Fig El the range of /ij is [0, 1 + x\. 
In the second case there is an extra restriction fa < 1 + 9 corresponding to the forbidden 
interval [1 + 9, 1 + x], and in the third case the restriction x — X < pi < 1 + 9, corresponds 
to two forbidden intervals [0, x — A] U [1 + 9, 1 + x]. 

For the upper cut hexagon the range of endpoints fi[ = m'Jk of lattice paths originated 
on the upper left side is [0, 1 + A + 6 — x] in all three cases of Figj5l while restrictions apply 
in the first case 9 — x < p,' { < 1 + 9, and in the second case /4 < 1 + 9. The upper and lower 
half correspond to an admissible rhombus tiling if and only if the endpoints ^ and a4 +1 _j 
coincide. It gives the identity 



which holds in all three cases. 

In the limit k — > oo the asymptotic of the counting formula of lattice paths gives the 
action functional 



where p e (p) = p{9 — p) and S\[p] is the action (12.31) . 

The limit density is the minimizer of this functional with the restrictions described above. 
The extrema of S[p] are solutions to the integral equation 



v'k+i-i = 1 + 9 - fa 



S[p] = S x [p] + S x +e-x[p} 



(3.1) 




subject to the restriction Jj p{u)du = 1. Here I = [Max(0, x — A), 1 + Min(x, 9)}. 
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Theorem 3.1. Depending on the geometry of the cut hexagon, the limit density has the 
following form. 

(1) The limit density p(z) is supported on the interval [a,b] C [Max(0,:r — A),l + 
Min(z,0)] with 

^X(X + 6-x) - y/xO(l + \ + 0) S ' 

\~+~9 

v / \(\ + 9-x) + y/x8(l + \ + 6) 



x + 9 

and on this interval is given by 



1 / , Ib — z ll + 9 — a . lb — z ll + x — a 

= - If 



(3.2) 



piz) = - Arctan W \- — - + Arctan 

Iy ' 7r \ V^-aV 1+0-6 \lz-a\l + x-b 

(3.3) — Arctan \ -\ - — Arctan \ \ 

\z-a\b V z-a\ X-x+b 

Outside of [a, b] the limit density is zero. 
(2) The limit density p(z) is 1 outside[a,b] i.e. on [Max(0, x — A), a] U [b, 1 + Min(x, 9)]], 
while inside [a, b] it is given by 



1/ A b-z 1+9-a . b-z 1+x-a 

piz) = H — Arctan \ \ — Arctan \ \ 

HK J 7T I \z-a\l + 9-b \z-a\l+x-b 



b — z fa . b — z X — x + a 

(3.4) + Arctan \ - Arctan \ \ - - 

v ' \z-a\b V z-a\l X-x + b 

(3) p(z) — 1, z G [Max(0, x — X),a], and p(z) = 0, z G [b, 1 + Min(x, 9)}} . In this case a 
and b are given by the same formulae as above while 



1 / , Ib — z ll + 9 — a . Ib — z ll + x — a 



piz) = - Arctan \ \ - — - + Arctan 

Iy ' 7r \ \lz-a\ll + 9-b \lz-a\l + x-b 

, n ^ x . Ib — z fa . Ib — z jX — x + a 

(3.5) + Arctan W w- — Arctan 



z — a\J b \ z — a \ X — x + b 

(4) p(z) = l,z e [Max(0,x - A), a], piz) = 0,z e [b,l + Mm(x,9)}. In this case 
l+ l ; , jb — z jl + 9 — a , , jb — z jl + x — a 



7T \ 



p(z) = H — Arctan \ \ — - — Arctan 



z — a\l + 9 — b \ z — a \ 1 + x — b 



b—z a b—z X—x+a 

(3.6) — Arctan \ * - — Arctan \ \ - 

" z-a\ b V z-a\ X-x + b 
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Proof. 1. Let us first assume that the density p is supported on an interval [a, b] with 
Max(0, x-A)<a<6<l + Min(x, 9), and that p(a) = p{b) = 0. We must solve eq. ([3"TTT) 

with the normalization p{u)du = 1. The corresponding resolvent reads: 
F{z) = vM™ ~ 1) f 1 dv 1 Log (ft + v){5 + v) 



2tt J y/ v (l- v ) w-v (7-v)(ti-v) 



i Log (p + w)(6 + w) + f / Arctan ]_ — w l_n_ + Arctan / 1 " 



2 (7 — w)(r] — w) \ \ w \ rj — 1 V u> y 7 — 1 



1-w p . 1-w 
—Arctan \ \ / — Arctan 



w V 1 + (3 ' V w \ 1 + 5 



where w = f— 2 and: 

b—a 



a 1+x — a \—x+a 1+9— a 

(3= l ' ^ = ~^ ' ° = ~^ ' ^ = ~u • 

b—a b—a b—a b—a 

The constrain on the density translates to the condition F(z) = 1/z + 0(1/ z 2 ), as z —>■ 00 
on the resolvent. Using identities ( 12.1 7ft ( 12T81) we obtain the equations: 



(v / ^TT+ y/p)(V6 + i + V6) = (v^+ v / 7^T)(v / ^+ v^F^) 

\VP + i - v^) 2 + Wi- Vi^) 2 + (Vs + T - V~5) 2 + (Vv- v^FT) 2 ) 



4 

Solving these equations for a and b we get the formulae (13.2ft . 

2. Next, let us assume the density is one near the boundaries, and varies only on an 
interval [a, b\. This happens, for example, for sufficiently small x, when < x < 9 (see the 
left figure on Figj5]). In this case p(z) = 1 for z E [0, a] U [b, 1 + x]. Eqn.f l3.ll) turns into: 

1 Q - a) 2 (l + x - p)(X - x + p) = f b p(u)du 
1 ' ] 2 ° g n{b-ii)*(l + e-ii) t p-u 

and we have 



= y/w(w-l) r 1 dv 1 ^ v 2 ( 7 - + v) 



27T 7 ^(1 - v) W - V ° (1-V) 2 (P + V)( V -V) 



- Log 7 W2( 7" w)(<y + W) , + i Arctan J J-^ . /3I - Arctan , / 1 " W ^ 



2 (1 — w) 2 ([3 + w)(t] — w) \ \ w y r] — 1 V w y 7 



1-w P . 1-w 5 

-Arctan \ / * / — Arctan 



u; V 1 + /^ V w V 1 + 5 
where 

a 1 + x — a A — x + a 1 + 9 — 

P = l ' 7 = ~^ ' 6= ~^ ' ^ = ~u 

b—a b—a b—a b—a 
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The constraint for the density implies F(z) = h ~ a ~ x + 0(-^ for large z. This gives: 



b-a-x = ^ U- + v^) 2 - (V7- V 7 ^!) 2 



Remarkably, this leads to the same result (I3.2p as above, and for the density we obtain 

dS2D 

3. In mixed cases where p(z) = 1 for z < a and p(z) = for z > b the resolvents are: 



_, ^w(w - 1) f 1 dv 1 v 2 (<f + v) 



27T 7 ^(1 - u) io - v (/3 + v) (7 - -u) (77 - i>) 



y/ W ( W -l) ^ dv 1 ( / 9 + V )( 7 _ + 



2vr 7 ^(1 - w ) u> - v (1 - v) 2 (r] - v) 

These formulae lead to the same formulae for a and b and to the formulae ( 13. 5ft (13.6j) for 
the densities. 

□ 

Remark 3.2. At the boundary intervals of length 1 (they correspond to cuts at x = and 
x = A + 9) the limit density is 1 . This can be derived from the formulae in the theorem and 
from the identities 



(3.8) 





- 1 = 


Max(xA(l + A + 6 


>),0(A + 0-x))* 


\fl - y/l 


- 1 = 


Min(xA(l + A + 9 


),6(\ + 6-x))* 




v 7 ^ = 


Max(A(A + 0-x) 


,x9(l + \ + 9))* 




= 


Min(A(A + fl-x). 


,x9(l+\ + 9))* 


Vv + Vv 


- 1 = 


Max(:rA,0(l + A - 


-9){\ + 9-x))* 


Vv- Vv 


- 1 = 


Mm(x\,9{l + A 4 


■6)(\ + 6-x))* 


Vs + i + 


Vs = 


Max(sc0,A(A + 0- 


-x)(l + \ + 9))* 


Vs + i- 


sT5 = 


Mm(x9,\(\ + 9- 


■x)(l + \ + 9))* 
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which, in turn, implies 



(3.9) 



i + e-a 



1 + 0-6 



1 + x — a 



1 + x — b 



X — x + a 



X — x + b 



xX + ^9(1 + X + 9){X + 9 



x 



X + 9 

xX- y/0(l + X + 9)(X + 9-x) 



X + 9 

^/xX(l + X + 9) + ^9{X + 9-x) 



X + 



y/xMX + X + 9) - ^9(X + 6-xj 



X + 9 



r x9- ^X{X + 9 -x)(\ + X + 




The general picture is summarized in FigEJ The density p varies only within the ellipse 
inscribed inside the hexagon. The points where the ellipse touching the boundary of the 
hexagon have x-coordinates x\ — and x 6 = X + 9 (vertical tangents), x 2 = j^j, x% = -^j, 

>w "' rH] anc l x 5 = X ^i+x^ ' wr th x\ < x 2 < x 3 < x 4 < x 5 < x e . They correspond, 
a, b = 1 + x, a = 0, b = 1 + 9, and a = x — X. 



x i 



1+0 

respectively, to b 

The limit density is constant outside of the ellipse Figj6] and is given by: 

(i) x G [xi, X2}: top and bottom p = 1 

(ii) x G [x 2 ,x 3 ]: top p = 0, bottom p = 1 
(hi) x G [#3, £4]: top and bottom p = 

(iv) x G [x^ £5]: top p = 1, bottom p = 

(v) x G [x 5 ,x 6 ]: top and bottom p = 1. 

The explicit expressions for p inside [a,b] are given by eq. (13.41) in the cases 
eq. fl3.5p in the case (ii), eq. fl3.3p in the case (iii), and eq.f l3.6p in the case (iv). 



and 



4. Comparison with the general results of Kenyon-Okounkov- Sheffield 
4.1. Equations for the slope. 

4.1.1. Limit shapes for rhombi tilings and more generally for dimer models on bipartite 
graphs are known to be solutions to certain variational problem [T|[5]. They minimize the 
rate functional which is determined by the free energy of the system on a torus. 

Assume we have a polygonal region with sides s\, . . . , s 3 d of lengths Ni, . . . , N 3d parallel 
to the axes of the triangular lattice, such that three possible axes of the lattice alternate 
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1 




(0 ; (M) ; (Mi) i ov) H 



P=i 






Figure 6. The full solution for the density p of path points (by vertical 
sections) in the case of an hexagon of size 1 x A x 9. The pink regions 
correspond to p — 1, while the blue ones stand for p = 0, both frozen. The 
domain of segments [a, b] over which p varies sweeps the inner ellipse. 

in the sequence Sj. Consider random tilings of such region with the probability oc q vol ( h ) 
where h is the height function corresponding to the tiling. When q — > 1 such that q~ ec 
and eiVj are kept finite random tilings develop the limit shape at the macroscopical scale 
[I]. Moreover in this case the limit shape is described by an algebraic curve which can be 
explicitly constructed from the geometric data of the region. 

Recall that rhombus tilings of a simply connected domain are in a bijection with the 
height functions (certain discrete surfaces in 3D over this domain). When the domain is 
not simply connected the height function does not exist but what would be its gradient still 
makes sense. It is called the slope of the tiling. The slope describes the local proportions 
of tiles of types a, b, c as follows. Consider the tiling of a toroidal domain of fixed size, let 
N a , Nb, N c be the numbers of tiles used in the tiling with N = N a + Nb + N c , then the slope 
is defined as (s,t) = Vh = (h x ,h y ) = (N a /N, N b /N) (here, x,y refer to coordinates in a 
basis (u,v), with u 2 = v 2 = 1 and u.v = —1). 

4.1.2. As explained in [5][1] there are two natural complex functions on the domain in 
question, z and w parameterizing the height function: 



(4.1) 
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In the limit when q — > 1 and the regions increase at the consistent rate (see above), these 
variables satisfy the system of equations: 

P{z,w) = 
Re(z x /z + Wy/w) = c 

for some constant c. In addition z and w should satisfy the boundary conditions determined 
by the geometry of the region and restrictions on tilings. 

The function P(z, w) is polynomial for periodically weighted lattices. It is the determi- 
nant of the Kasteleyn matrix with quasiperiodi<£| boundary conditions which is the partition 
function of the dimer model on the (toroidal) fundamental domain of its weights. 

In our case 

P(z, w) = 1 + z + w. 

Solving the system for z and w by complex characteristics we change it to the system of 
equations 

(4.2) P(z,w) = 

(4.3) Q{e- cx z,e~ cy w) = 

The function Q is an analytic function, which depends on the boundary of the region. 

When the region is polygonal, the function Q is polynomial. It can be computed accord- 
ing to the description of [I] (Section 3.3.3). 

One of the remarkable phenomena we observe in tilings of large domains is the existence 
of the frozen regions near the boundary. In these regions the height function is linear. 
In the case of rhombus tilings of a hexagon, the curve separating the frozen part and the 
disordered phase (where the height function varies smoothly) is an ellipse [2], and it was 
called the arctic circle. In [4j such curves are called clout curves. One of the important 
corollaries is that for polygonal regions such curves are real algebraic in e cx ,e cy . It is also 
important that they form a family of curves as one varies c. 

When c — > 0, the second equation in (14. 2p becomes Q(z, xz + yw) = 0. 

The arctic circle (the boundary of the limit shape) in the case of polygonal region is 
an algebraic curve intersecting each linear piece of the boundary once where it is tangent 
to the boundary. It may develop cusps. It may also be disconnected, i.e. it may have 
"islands" and bubbles" (in the terminology of |4j), see Fig. 5, fig. 17 from [4]. 

4.2. Case of the cut hexagon. As explained before, the limiting density of exit points of 
LGV paths on a large cut hexagon kxnxn with n/k = A fixed, k — > oo is identical to that 
of passing points of LGV paths through the main diagonal (which is where we make the 
cut). These paths describe the rhombus tilings of the hexagon which is obtained from the 
cut hexagon by gluing to the reflection w.r.t. the cut. In the language of [I], this density 

1 Quasiperidicity means there are extra factors z^ 1 and u> ±:L for the dimers crossing two given lines 
parallel to the x and y directions respectively, depending on the Kasteleyn orientation of the corresponding 
edges. 
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is the density of tilted rhombi and is determined by the limit shape of rhombus tilings of 
the complete hexagon. 

For the unweighted tilings by rhombi we have P(z,w) = 1 + z + w. The polynomial Q 
for such tiling of the hexagon has degree 2 ( half of the number of sites of the region) and 
reads: 

Q(z, xz + yw) = X 2 (z + z 2 ) — (xz + yw) 2 — (1 + X)(xz + yw) 
The equations (14.21) for c = become quadratic equation for z: 

(A 2 - (x - yf)z 2 + (A 2 - (1 + A - 2y)(x - y))z + y(l + A - y) = 

The arctic circle is the curve in (x, y) on which the discriminant of this equation vanishes: 

= (A 2 -(l + A-2 2/ )(x- 2/ )) 2 -%(l + A- Z/ )(A 2 -(x-y) 2 ) 

= (1 + 2\){x - yf + A 2 (a; + y - A - l) 2 - A 2 (l + 2A) 

It is tangent to the six edges of the hexagon, namely y = 0,y = X + l, x = 0,x = X + l, and 
y = x±X. This curve intersects the reflection axis y = x (the cut) at points x = A±l-|-yIi±2A 
which coincide with the ends of the support of p . 

The passing points of the LGV paths along the cut correspond to tiles of type a or b. 
The local density of such tiles is therefore (N a + N b )/N = h x + h y . Hence the density (13.31) 
should be compared to the "vertical" slope h x + h y in our solution. 

Solving for z at y = x, we find that A 2 (^ + ±) 2 = (x - ^) 2 - so that: 



-1 ± 



(x — a)(x — b) 



ab 



'x — a)(x — b) 



ab 



1 . (b — x)(x — a) 

h x = h v = -Arctan ■ 



and finally we get 

h - I 

IL x i»y ii , 

7i V ab 

and therefore p(x) = h x + h y as expected. 

4.3. Case of the full hexagon. Here we will compare the results of previous sections 
with the computation of the limit shape based on the equation ( 14.21) . 

In the case of the full hexagon with edges 1, A, 9 the equations (14.21) becomes 

= 1 + z + w 

= \6(z + z 2 ) - {xz + yw) 2 - (A - 9)z(xz + yw) - (1 + X)(xz + yw) 

The arctic curve appear as the condition that the discriminant of this equation vanishes. 
It is the ellipse: 

((1 + A)x + (1 + 8)y - (1 + A)(l + 9)Y ((1 + X)x - (1 + 8) V f 

(«) A TTaT^ " + * Ye " + W + "» 
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tangent to respectively x = 0,x = 9,y = 0, y = X, y = x + X and y = x — 9. 

To get the density of LGV points, let us solve for z on some arbitrary vertical line 
y = x + t. It is easy to see that the line crosses the arctic curve at the points 

\/\(t + 9)T y/6(\-t)(l + \ X ' 



x± 



x + e 




These are identified with the ends of the support of p given in eq. fl3.2p . namely x_ = a and 
x + = b, upon the substitution x = \ — t. The solution for z reads: 

z = ~ 2 (t + 9)(X-t) ( t +( X - t - x )( t + e ) + ( X - t )( t + x ) ± ( X + °) \/{x-a){x-b)) 
and we finally get 



h - Cretan. ( X + OVW^)^ 



7T (X-t-x)(t + 9) + t + (A -t)(t + x) 



1 (X + 9)^(b-x)(x-a) 

h„ = — Arctan ■ 



it (X-t + x)(t + 9)-t-(X-t)(t + x) 



Using the addition formula for Arctan, we finally get 
(4.5) h x + h y = ±Arctan (A + &) ^ ~ ^ 



7T \e - t{l + 9) + 2x(x + t-i-Nf) 

To compare this with the results on the LGV point densities, let us assume for definiteness 
that we have x 3 < x < x± (case (iii)), namely — 9<t< A^. Using eqs. (13.21) and (13.91) 
and the restriction on t, we find that 



11 + 9 


— a 


1 + 6 


- b 




fa 




\jb 


'l + x 


— a 


1 + x 


-b 


X — x + a 


X — x 


+ b 



^9(9- 


f t)(l + A + 


) + VK X - 1) 


y/9(9- 


ft)(l + A + 


) - a/A(A - t) 


V 9(X 


-t)(l + X + 9 


)- y /X(9 + t) 




-t)(l+X + 6 


>) + ^9(9 + t) 


Va(a 


-t)(l + X + 6 


)) + ^9{9 + t) 


y/X(X 


-t)(l + X + 6 


>)-^9(9 + t) 


y/X(9 + t)(l + X + 9 


) - >/e{\ - t) 


^X(9- 


+ t)(l + X + 9 


) + V0( x - 1) 



We use the subtraction/addition formulae for Arctan to get: 

y/(b-z)(z-a)(X + 9) 



Arctan \ I \ I — — ; — Arctan \ I \ - — Arctan 







jl + x- 


a 




a \ 


' 1 + x- 


b 


'b- 


z 


jl + 9- 


a 



b 9X-t(l + 9) + z(X-9) 



Arctan t/- -»/- - A retain/- — X + a = Arctan V ' ''' U ' 



z~a\ 1 + 9-b V z-a\ X-x + b 9X + t(l + 9) — z(X — 9) 
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Figure 7. A cut hexagon with zig-zag fixed boundary. 



and finally 



p(z) = 1 Arctan f> ~ + g) + Arctall . V<" - ~ ^ + 



7T v ex-t{\ + e) + z(\ -e) ox + t(i + 0) - z{x - 0) 

1 v /(6-^)( 2; -a)(A + 

— Arctan ■ 



7T ^A-t(l + ^) + 2^(z + t- 1 



A+6 
2 



in agreement with (14. 5p . 

The other cases follow analogously. 

5. Symmetries 

5.1. Case of a cut hexagon with a fixed zig-zag boundary. 

5.1.1. The "half cut hexagon (see FigJTJ), of the size k x 2n x In + k is another case when 
the total number of lattice path configurations nicely factorizes. In this case the lattice 
paths start at each point of the left boundary and ends on the right boundary. The low 
boundary is a horizontal zig-zag line. 

The LGV formula expresses the total number of configurations of such non-intersecting 
paths, which start at (0, 2i — 2), i = 1,2, ...,k, and end at (2n, 2rrij — 2), j = 1,2, ...,k, 
1 < rrij < k + n as the determinant: deti<jj<jfc(p^.). Here p\ n ^ is the numbers of paths 
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from (0, 2% — 2) to (2n, 2m — 2). By the standard reflection principle, this number pf^ is 
exactly the number of unrestricted paths from (0, 2i — 2) to (2n, 2m — 2) minus the number 
of of unrestricted paths from (0, —2i) to (2n, 2m — 2): 

2n \ ( 2n 
,n + m — ij \n + m + i — 1 



in) 
Vim 



Finally we have the factorization formula 



A (2n + 2*-2)! 



(5.1) Z fc>B ({m})= det = (rfX)=II 



i<< J<fc ' mj 1 4 (" + wii + k - 1) ! (n - m, L + fc) ! 

1 = 1 

1 I ( m i _ m ^> 1 i ( m * + m j _ !) 

l<i<j<k l<i<j<k 

Let us now compute the asymptotic of this number, when k — > oo, and say n/k = a 
fixed, while p^ = TOj/A; converges to a continuous function on [0,1], or equivalently gets 
distributed on [0, a + 1] with the limit density p(p). 

We have 

1 1 k 

(5.2) - -Log(Z fe , n ({m})) = C + - ^((a + 1 + in) Log(a + 1 + &)+ 

i=i 



(a + 1 - pi) Log(a + 1 - ^)) - A LoglpJ - P? 



A; 2 

for some constant C depending only on a. As k — > oo this gives the action (the rate 
functional) 

(5.3) 

ra+l 

S[p] = C + / ((a + 1 + p) Log(a + 1 + p) + (a + 1 - p) Log(a + 1 - p))p(p)o?p- 

/■a+l /-a+l 

/ / Log|p 2 - v 2 \p{n)p{v)dndv 
Jo Jo 

Fixed points of this functional are solutions to the integral equation 

a. + 1 + p r +1 / i i \ 

Log — — = 4 + — — a(u)du 

a + l — p Jo \p — w p + w / 

subject to the constraint J a(u)du = 1. As before assume the density is zero for z e 
[a, a + l]. It is convenient to extent the problem to the symmetric interval [—a, a] assuming 
that a(z) is extended as a symmetric function (7(2) = a{—z). The integral equation becomes 



a + l +fl /•« l 

Log — — = -f cr(w)dw 

a + 1 - p J_ a fi-u 
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with the constraint J^ a a(u)du = 2. 

Similarly to previous sections it is easy to obtain an explicit formula for the resolvent 

= Il a ^{u)du: 



7r J a ^ u 2 - a 2 a + l- uz-u 

The constraint on the density is equivalent to the condition G(z) — | + as z — > oo. 

Solving this, we arrive to a = a + 1 — ^ 3 + 4q . Combining this with the explicit evaluation 
of the integral for G(z) we obtain the limit density: 

,r as t x 2 a V3 + 4a - 4^ 2 

5.4 erm = — Arctan 

V ; V ; 7T l + 2a 



for z G [0, y | + a] and cr(;z) = for z > yj\ + a. 

5.1.2. We could have found this result a priori, by recognizing that the half-cut hexagon 
from the previous section is the fundamental domain for the obvious reflection symmetry 
of the cut hexagon from the first section. 

Because the limit density is the unique minimizer of the rate functional, it is invariant 
with respect to this reflection. In terms of the limit density or the cut hexagon from section 
CD this is the symmetry p(z) = p(X + 1 — z). 

Taking into account the change of variables from Fig. [I] to Fig. [7] we can see that indeed 

p(z + a + 1) = a(z) 
where A + 1 = 2(a + 1), p(t) is given by (11.151) and a(z) is by (15.41) . 
5.2. The reflection principle for limit densities. 

5.2.1. The reflection principle which was illustrated above can be used to compute the 
limit density for a random tiling along the free boundary in a more general setting. In 
particular, it can be used to compute such density from the limit shape on doubled domain 
with fixed boundary. 

Assume that the free boundary of a domain D is a collection of segments on a line and 
let r D r be the reflection with respect to this line. The doubled domain D U T> r is closed 
with feasible (in a sense of [1]) boundary. It is invariant with respect to the reflection 
D — > D r . If we extend the definition of weights of tiling to D r such that the weights will 
also be invariant with respect to this reflection, the resulting limit shape density will also 
be invariant. It is also clear that in this case the limit shape density for tilings of the 
domain D U D r at the symmetry axis is the same as the the limit density for the tilings of 
T) with the free boundary condition along this axis. Thus, the problem of computing the 
density of tiles at the free boundary reduces to the computation of the limit shape for a 
symmetric domain with fixed boundary along the symmetry axis. 

As an example we can compare the limit density at the diagonal of the symmetric hexagon 
with A = 9 with the density at the free boundary of the cut hexagon with the same value of 
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(a) (b) (c) 

Figure 8. The regular hexagon is decomposed into six triangles (a). The 
two inner boundaries are both parts of lines that divide the hexagon in two 
(dashed lines), each corresponding to cuts along lattice paths made of the 
two indicated tiles, for which the density along the cut is known. In (b) 
we interpret the tiling of the triangle in terms of lattice paths starting at 
the fixed edge and ending on the two free inner boundaries, along the third 
type of edge, not used for the paths in (a). In (c), the same tiling is now 
interpreted differently in terms of paths that start on the lower free boundary 
via tiles of the third type, and end on the top one with the two first type of 
tiles. 

A. The symmetry axis of the hexagon correspond x = A in the computation of the density 
along slices of the hexagon. The density along this line is given by eq. fl3.3p with x = 9 = A, 
and with a and b given by eq. fl3.2p . These expressions for the hexagon case are identical to 
formulae fll.l3pfll.l5p for the density of endpoints of lattice paths along the free boundary 
for the half- hexagon of the same shape. 

5.2.2. Now let us discuss limit shapes corresponding to lattice paths in the cut hexagon 
with forbidden intervals along the cut. 

Lattice paths and equivalent to tilings by rhombi which is also equivalent to monotonic 
piles of unit cubes in the corner of a "3D-room" . This is simply another way to phrase that 
rhombus tiling configurations are equivalent to the corresponding height functions. 

The case of the cut hexagon with one forbidden interval corresponds to pile of cubes 
in a rectangular room with two parallelepipeds removed from it (see FigJHl (a)). One is 
removed from the upper closest corner, one is removed from the lower corner of the other 
side. They should have the complementary heights (their sum is equal to the height of the 
room) and they projections on the (x,y) plane should intersect. If this configuration is 
symmetric with respect to the reflection (x,y) — > (y,x), the lattice paths corresponding to 
such a pile of cubes will have the forbidden interval along the diagonal of the square which 
is the intersection of the lower side of the upper parallelepiped and the upper side of the 
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Figure 9. The 3D realization of forbidden (a) and fully packed (b) bound- 
aries. Both are obtained by removing from the large parallelepiped (the room, 
in dashed lines) two smaller ones touching opposite corners of the room. In 
the case (a), they share a plane domain (a square in the case of regular ob- 
jects), while in case (b) they intersect along a segment. The diagonal of the 
square and the segment, in projection, correspond to respectively the forbid- 
den interval and the fully packed one. We have also represented the expected 
limit shapes for the associated rhombus tilings. 

lower parallelepiped. If we cut this rectangular room along the diagonal, the corresponding 
lattice paths will be exactly lattice paths in the cut hexagon with a forbidden interval along 
the cut. 

In this case the density of the LGV paths along the cut is determined by the limit slope 
of the corresponding tilings, or equivalently by the limit height of corresponding piles of 
cubes. The boundary of the limit shape in this case is disconnected. In terminology [4J 
this limit shape is the result of the degeneration of the island with the hexagon inside and 
six cusps. The degeneration "flattens" the hexagon to a rhombus with horizontal main 
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diagonal (two vertical sides shrink to points). As a result six cusps around the hexagonal 
island degenerate into four cusps around the rhombus. 

A similar interpretation holds for the density of end points of lattice paths with fully 
packed interval along the free boundary of the cut hexagon. In this case piles of cubes 
also should be in the a 3D rectangular room with two parallelepipeds removed from it (see 
Figj9] (b)). As in the case of forbidden intervals, one of them is removed from the upper 
closest corner, one is removed from the lower corner of the other side. The difference is 
that their projections on the (x, y) plane should be complementary (they should intersect 
over exactly one point). The sum of heights should be greater then the height of the room. 
In this case parallelepipeds intersect over a vertical segment. When such a configuration 
is symmetric with respect to the reflection (x,y) — > (y,x), the piles of cubes in such room 
with removed parallelepipeds are equivalent to the lattice paths on the left side continuing 
to the lattice paths on the right side with the fully packed interval at the diagonal cut (the 
interval where two parallelepipeds intersect). 

In this case the limit shape will be discontinuous. In the terminology of [I] this case is 
also a degeneration of an island but the hexagon now becomes "squashed" sidewise such 
that it degenerates to a vertical segment (tilted sides degenerate to points). In this case 
six cusps also degenerate to four cusps. 

5.2.3. This symmetry principle extends to situations with other symmetries. Consider for 
instance a regular hexagon (corresponding to A = 9 = 1 in the previous section). Then the 
domain is invariant under rotations by 27r/3, and reflections w.r.t. axes joining opposite 
vertices. Let us cut the hexagon into six triangles according to these symmetries, and focus 
on one of them (see Fig 18] (a)). 

The counter clock- wise rotation by 27r/3 brings a tiling to a tiling. It maps the square 
looking tiles from the upper side of the triangle from Fig. 8(c) to the tiles along the lower 
side of the triangle from Fig. 8(b) where lattice paths end. 

Because the hexagon is invariant with respect to rotations, because the action function 
(the entropy, or the rate function) is invariant with respect to such rotations, and because 
it has unique minimizer the limit hight function is also invariant with respect to such 
rotations, if one takes into account that the rotations also rotate tiles. 

In particular , if a(z) is the density of lattice paths along the lower boundary of the 
triangle from Fig. 8(c) and p(z) is the density of paths along the upper boundary of the 
same triangle, we have p(z) = 1 — cr(z). Here z is the distance from the right lower corner. 

The density p(z) we computed earlier in eq. (11.151) for all values of A. For A = 1 we have: 



(5.5) 




ze[0,i-f] 

26[l-f,l] 



Here z is the coordinate along the upper side of the triangle with z = being the upper 
corner and z = 1 being the lower right corner. 
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For the density a(z) we obtain 




z e [o, l - & ] 



( 5 - 6 ) a ( Z ) = S 1 Arrt.an I 1 , <= U - & 



, ZE fl-^,11 

More generally, from the results for lattice paths in a hexagon, we obtain the density of 
the lattice paths in the triangle from Fig. 8(c). Let L x be the line parallel to the upper 
side of the triangle which is at the distance x from the left lower corner. In the formulae 
below we use the coordinate z which is at the left side of the triangle and changes from 
to x along L x as we move downward. Let p{z,x) be the density of paths crossing L x , 
0<z<x. 

For 0<x<l — ^we have p(z, x) = 1 for all < z < x. 
For 1 — ^ < x < | and z E [0, a] we have p(z, x) = 1 and 




. . , .. . .1 + 9 — a . /b-z jl + x 

p(z, x) = 1 H — Arctan \ \ — Arctan 



9 — b \ z — a \ I + x — b 



b—z a b—zX—x+a 

(5.7) + Arctan a / W- — Arctan \ \ - 

Xl z-aV b V z-oV X-x + b 



when z G [a, x]. Here a = ^ - |a/(2 - x)3a; and 6 = + 1^(2 - x)3x. 
For | < x < 1, we have p(z, x) = when < z < a and 



p(z,x) = — I Arctan W \l — — -+ Arctan 




z — a \ 1 + x 



, ., „la . b — z IX — x + a 

(5.8) — Arctan \ * - — Arctan \ \ 

" - b V z - a\ X - x + b 

These formulae compare with eq. (15. 5p as p(z) = p(z, 1). 

5.2.4. Limit shape of TSSCPP. The Totally Symmetric Self-Complementary Plane Parti- 
tions (TSSCPP) are in bijection with the tilings of a half equilateral triangle of edge size 
2n, with fixed straight boundary along one edge, fized zig-zag boundary along the median 
perpendicular to that edge, and free boundary along the remaining edge (see Fig jiu]) . 

The region in this problem can be naturally identified with the fundamental domain of 
the full, fully symmetric hexagon. It is natural to choose it as the triangle with < x < 1 
and x < z < ^rp. Under this identification the zig-gag line becomes the line z = x. As 
n — ► oo it becomes the symmetry axis of the limit shape of the height function for the 
equilateral hexagon. 

Using the same symmetry arguments as above we can compute the limit density of paths 
crossing the vertical line L x which is at the distance x from the left side of the triangle. 
Denote this density by a(z, x). 
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Figure 10. TSSCPP are the rhombus tilings of a half-equilateral triangle, 
with fixed boundary along the half-edge and the zig-zag median, and free 
boundary along the remaining edge. We have indicated with dots the starting 
and ending points of the non-intersecting paths with a floor, in bijection with 
the tilings. The exit tiles are those marked by dots along the free boundary. 
We expect the zig-zag line to become the reflection symmetry axis of the full 
equilateral triangle in the continuum. 




When 0<x<l — ^we have a(z, x) = 1 for < z < a and a(z, x) is given by eq. (15.71) 
for a < z < 

When 1 — ^ < x < 1/2 the density is given by eq. (15. 7p for all x < z < . It is given 
by eq. (15 .7p for all x < z < when 1/2 < x < 1. 

These expressions agree with the function eq. (15. 6p : a(x, x) = a(x). 
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5.3. Some other examples. Another example of the application of the symmetry prin- 
ciple is the computation of the limit shape for the equilateral hexagonal region with free 
boundary condition on two adjacent sides. In this case one should glue three such hexagons 
into the region shown on FigJTTTa). The limit shape for this region can be found using the 
methods of [I]. 

The cut hexagon with both sides having free boundary conditions on tilings correspond 
to a random tiling of the infinite region obtained by gluing infinitely many hexagons along 
vertical boundaries FigJTTTb). This case also describes the hexagon where two opposite 
sides have the free boundary. 

For more complicated examples of regions with free boundaries the corresponding re- 
gion with fixed boundary may not be simply connected (see Fig JTTT c)). nor planar (see 
FigjTT](d)). 
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